Freezing of quantum correlations under nondissipative decoherence is universal 



o 

(N 

CI. 
< 
CO 



Oh 



> 



o 

CO 



X 



Benjamin Aaronson/ Rosario Lo Franco,^ and Gerardo Adesso^ 

^School of Mathematical Sciences, The University of Nottingham, 

University Park, Nottingham NG7 2RD, United Kingdom 

^ Dipartimento di Fisica e Chimica, Universitd di Palermo, via Archirafi 36, 90123 Palermo, Italy 

(Dated: April 3, 2013) 

We show that the phenomenon of frozen discord, exhibited by specific classes of two-qubit states under local 
nondissipative decoherent evolutions, is a common feature of all known bona fide measures of general quantum 
correlations. All those measures, despite inducing generically inequivalent orderings on the set of nonclassically 
correlated states, return a constant value in the considered settings. Every communication protocol which relies 
on quantum correlations as resource will run with a performance completely unaffected by noise in the specified 
dynamical conditions. We provide a geometric interpretation of this phenomenon. 
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Quantum correlations, seminally quantified by the quantum 
discord OEl. stand as one of the most general manifestations 
of nonclassicality in composite systems. They can be revealed 
in the process of locally measuring a subsystem, even in states 
where entanglement or nonlocality are absent. Despite a mas- 
sive surge in recent studies investigating interpretation, quan- 
tification, and applications of discord and related quantifiers 
of quantum correlations |3 1, it is a fact that these quantities 
remain far less understood than entanglement fT] . Few proper- 
ties are now set in stone as necessary requirements to identify 
a quantifier, say (3, as a bona fide measure of general quantum 
correlations, revealed on the subsystem A, in bipartite states 
Pab- We list them below: 

PL Vanishing on classical-quantum (CQ) states: Q(pab) > 
for all Pab and Q(xab) = for all states xab belonging to 
the set QcQ of CQ states, xab = Zz/^/l^'X^'U ^ t/^, with {/?,} 
being a probability distribution, |/)a an orthonormal basis for 
subsystem A, and T/g arbitrary states for subsystem B; 

P2. Invariance under local unitaries: Q{(Ua Ub)Pab(Ua ^ 

Ub)^) = Q{PAB)\ 

P3. Nonincreasing under local operations on the unmea- 
sured party B\ Q{(1a <8) ^b)[Pab]) < Q(Pab), where Og is a 
completely positive and trace preserving map (i.e., a quantum 
channel) acting on subsystem B; 

P4. Reduction to entanglement on pure states: 
QQ^abX^abD is an entanglement monotone IH. 

These criteria are certainly not complete, as in particular a 
definition of the operations on A, or of the allowed classical 
communication between A and B, under which discord and 
any valid Q should not increase, has not been accomplished 
ElllHIl- In the quest to unveil the most essential signatures 
of quantumness in complex systems, it is then of wide inter- 
est to identify physically insightful properties that underly the 
notion of quantum correlations as opposed to entanglement, 
and should then be reflected by any valid measure thereof. 
One such property is, for instance, the absence of monogamy 
I181I . This area of investigation has also a strong technologi- 
cal motivation |9|: it is believed, and in some case proven, 
that states with quantum correlations other than entanglement 
can be employed as resources for several quantum computa- 



tion |[TQV[T21 . communication lfT3HT6ll , and metrology setups 
(TTIITHI. Identifying the distinctive traits of such correlations, 
in particular with respect to their resilience under noise, might 
lead to valuable recipes for their practical exploitation. 

Numerous works have in fact investigated the dynamics of 
general quantum correlations in open quantum systems un- 
dergoing various types of Markovian or non-Markovian evo- 
lutions, as recently reviewed in O [191 1201 • One evident fea- 
ture is that discord is typically more robust than entanglement 
and does not suffer from sudden death issues Ii21|j22|; dis- 
cord can even be created by local operations on the measured 
party A (nonunital channels if A is a qubit) ||6||23l. A partic- 
ularly fascinating phenomenon can occur for two-qubit states 
undergoing nondissipative decoherence: their discord can re- 
main constant, or frozen, for an interval of time in Markovian 
conditions |24|. A forever frozen discord 1.25. L or multiple in- 
tervals of recurring frozen discord |[26l - l28l . can further occur 
when the dynamics is non-Markovian. Necessary and suffi- 
cient conditions for the freezing have been derived in |29|. 
However, to date such a feature seemed bound to the choice 
of particular (mainly entropic) quantifiers of discord |3|. It is 
natural to question whether freezing just happens as a mathe- 
matical accident, or whether it bears a deeper physical mean- 
ing which should manifest independently of the adopted mea- 
sure. Answering this question is the purpose of our work. 

We consider a selection of essentially all the known bona 
fide measures of general quantum correlations defined in re- 
cent literature, and we find that they all freeze under the same 
dynamical conditions. A geometric analysis, inspired by 1301 , 
is carried out in order to provide a satisfactory interpretation to 
the universality of this dynamical phenomenon. Our conclu- 
sions demonstrate that for all quantum information protocols 
(e.g. Cni [IMSI im [m EIHSI) whose performance refies 
on some form of discord between two qubits, however quanti- 
fied, there exist unique noisy evolutions in the state space for 
which coherence and thus quantumness in the correlations are 
exactly preserved. The corresponding protocols will then run 
with efficiency unaffected by such noisy conditions. 

We focus our attention to Bell diagonal (BD) states of two 
qubits, i.e., states with maximally mixed marginals Bll34ll. 



Their density matrix can be written in Bloch form as p^^ = 
\(^AB + ZL ^i^ ^ ^f ) where {crf'^} denote the PauH ma- 
trices and the vector c = (ci,C2,C3) completely specifies the 
state. As the name suggests, BD states p^^ have the four 
maximally entangled Bell states as eigenvectors, with eigen- 
values A^^ = \[1 -h (-Ifci - (-1Y^^C2 + (-1)^C3], where 
a,b = 0,1. The conditions A^^^ > impose constraints 
on the entries of c, so that physically allowed BD states 
can be represented as points within a tetrahedron of vertices 
(-1, -1, -1), (-1, 1, 1), (1, -1, 1) and (1, 1, -1) in the three- 
dimensional space spanned by (ci, C2, C3) [34l. This geometric 
picture is very appealing to visualize dynamical trajectories 
1301 . For later convenience, we define an auxiliary ordered 
vector ^ = (^1 = |c/J,^2 = \ci2lg3 = \ci,\), where {/i,/2,/3} is 
a suitable permutation of {1, 2, 3} such that ^1 > ^2 > ^3- BD 
states with zero discord correspond then to ^2 = ^3 = 0- 

We study one-way measures Q of discord, that is, those 
which reveal quantumness of correlations as perceived by an 
observer probing only qubit A. Incidentally, whenever sym- 
metrized versions of these measures are available |3 1, they co- 
incide with their one-way counterparts for BD states, thus en- 
larging the scope of our study. The adopted measures, which 
include those reviewed in |3 1, others taken from recent litera- 
ture, and some defined or calculated here for the first time, are 
listed in the following [35 1. Most of these measures do not 
enjoy a closed analytical form in general, but are computable 
for BD states. Notice that pairs of distinct quantifiers selected 
from our list will generically induce inequivalent orderings on 
the set of non-CQ states, even within the BD class |36|. 

Quantum discord D\l,2\. Discord can be defined as the min- 
imum diff'erence in total correlations between the state p and 
the CQ state obtained after an optimized projective measure- 
ment {Uf} on A, D(pab) = min^nf } [KpAByKI^i nfpA^nf )]. 
Here I(pab) = S(pa) -^ S(pb) - S(pab) is the mutual infor- 
mation (371, with S(p) = -Tr(plog2p) being the von Neu- 
mann entropy. For BD states, S(p^) = S(p^^) = 1 while 
^(Plfi) = - 2a,Z;=o ^L log2 W^)- The discord can be com- 
puted as (301IS1 ^(pL) = 1 - ^(pL) + i^(^), where 
H(s) = -slog2 s-(l- ^)log2(l - s). 

One-way quantum deficit A^ |3^|40l. This measure quanti- 
fies, in a thermodynamical framework, the minimum entropy 
production after a projective measurement on A, A~^(pab) = 
min|nA} [S(Yji^fpAB^f) - S(pab)]- It can be interpreted as 
the amount of information in the state pab which cannot be 
localized via a one-way channel of classical communication 
from A to B. For BD states, A~^(Pab) = ^(P^b"^ I38|. 
Relative entropy of discord Dr |40-42J. This quantifier cap- 
tures the distance, as measured by the quantum relative en- 
tropy 1431, between pab and the set of CQ states, Dr(pab) = 
^^XABeacQ^(pAB\\XAB), where S(p\\cr) = Tr[p(log2P - 
log2 cr)]. For all bipartite states, Dr(pab) = ^~^(Pab) 1401 l42l . 
Adjusted geometric discord Dq. The geometric discord Do 
was defined in ||44l |45l as the minimum squared Hilbert- 
Schmidt distance from the set Qcq- However, it is now ac- 
knowledged that such a definition is flawed I.46.L as Dg can be 



modified arbitrarily by reversible operations on B. Moreover, 
Dg does not respect P3 |47, 48 1. It is customary, when dealing 
with the Hilbert-Schmidt distance, to normalize it by the pu- 
rity of the state, in order to correct the metric for the eff'ective 
dimension of the Hilbert space ||49| |50 1 . We then define the 
adjusted geometric discord as Dg(pab) = 2(min;^^^eacQ WPab- 
;rAfill2)/Tr(p^g), where \\M\\l = Tr[M^M]. Properties and ap- 
plications of this measure are discussed in [51 1. For BD states, 

using m we have DGipi^) = ^(^ + ^l)/(^ + ^1 + ^2 + ^3)- 
Trace-distance discord Dx |[52j[53l. As the name suggests, 
this measure is given by the minimum trace distance from the 
set of CQ st ates, D i(pab) = min^^^^eQcQ WPab - XabWu with 
||M||i = Tr[ Vm^M]. For BD states, such a measure is com- 
putable |[52j|54l[55| and one has simply Di(p^^) = ^2- 
Negativity of quantumness Qn 1 52 , 56 1 This measure corre- 
sponds to the minimum negativity (an entanglement measure 
1 57 1) created between the system AB and an apparatus C dur- 
ing a local projective measurement on A, according to the for- 
maHsm of |[56l|58l|59l. Referring the reader to (521 El for de- 
tails, we recall that if subsystem A is a qubit (as in our case), 
then Qn(Pab) = Di(pab)- 

Bures-distance discord Db. We can consider the Bures dis- 
tance from the set of CQ states as another measure of dis- 
cord, in analogy with the Bures measure of entanglement 

EQl. We then have Db(pab) = [(2 + V2)(l - V^max(PA5))] ' , 
where F^^^Qjab) = max^^^^^acQ F(pab,Xab) and F(p, cr) = 
{Tr[{ ^cr^/py]f is the Uhlmann fidelity f6ll. For BD 
states, in this paper we obtain an analytical expression for 
^max ll62l (see 16311 for an independent yet related deriva- 
tion). We notice that, for a given p^^, the CQ state which 
maximizes the fidelity is in general diff'erent from the one 
which minimizes trace distance, Hilbert-Schmidt distance, 
and relative entropy. Explicitly, we have i^maxCp^g) = ^ + 

\ msLX^ij^k) [ V(l + ^^■)^ - (^j - ^kf + ^J{l-Cif-\cj + Cky'\ 
where (/, j, k) denotes cyclic permutations of {1, 2, 3}. 
Fidelity -based measure Dp |6|. As a simple rescaling of 
the previous quantity, we also pick the geometric measure of 
quantumness discussed in (6115811641, which in our notation 
reads DpipAB) = 2[1 - i^max (Pas)], and is computable for BD 
states using the expression for Fmax given above. 
Local quantum uncertainty *tl ifTSl . The last measure we 
adopt, which is also the most recently introduced, quanti- 
fies the minimum quantum uncertainty on a single local ob- 
servable, IKpab) = ^^Ka Iwy(Pab, Ka 1b), where Ka is 
a Hermitian operator with nondegenerate spectrum acting on 
subsystem A, and /wy(p, K) = TrlpK^ - ^/pK ^/pK] is the 
Wigner-Yanase skew information |65|. The quantity tl is 
computable for all two-qubit states 1 18 1. Interestingly, we find 
that for BD states t((p^^) = Df(p^^), although the two mea- 
sures do not coincide in general. 

We consider local Markovian nondissipative decoherence 
channels acting independently on each of two qubits m = A,B 
in the state pab- The single-qubit Lindblad operator has the 
form £klPm] = ^{cr^PmO'^ - Pm), where y is the decoherence 
rate and k = 1,2, 3 denote bit flip, bit-phase flip, and phase 
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FIG. 1: (Color online). Time evolution of various measures of quantum correlations for BD states with initial condition c^^^(O) given by Eq. |2]), 
subject to local bit flip (k = 1), bit-phase flip (k = 2) or phase flip (k = 3) channels with rate y. In the range of parameters such that 2yt < - In c, 
depicted as a shaded (yellow online) region on the horizontal planes in each panel, all the considered measures remain frozen to their respective 
initial values. The plotted quantities_^are: (a) quantum discord D, equal to one-way quantum deficit A^, equal to relative entropy of discord 
Dr, (b) adjusted geometric discord Dg', (c) trace-distance discord Di, equal to negativity of quantumness Q^', (d) Bures-distance discord Db', 
(e) fidelity-based measure Dp, equal to local quantum uncertainty "U. The analytical expressions are given by combining Eq. ^ with Table |l| 



flip (alias phase damping) channels, respectively ||2T]|24I. We 
consider BD states p^y as inputs. The time-evolved states 
under the considered local channels remain in BD form, with 



Ci,Mt) = Cij(0)e-'^\ Ck(t) = Ck(0) 



lyt 



(1) 



where k selects the channel as explained above. Diff'erent 
initial conditions c(0) lead to varied dynamics of nonclassi- 
cal correlations, without any general agreement between the 
measures we consider. However, for each selected channel, 
indexed by k, we can choose a specific subset of initial condi- 
tions c^^^(O) 1 24 1, given by 

cfifd) = ±1, cfifd) = =Fcf (0), with |cf (0)1 = c (2) 

(notice that taking into account the sign freedom and the 
permutation of /, j i^ k there are four possible choices per 
channel). These conditions, which are equivalent to impos- 
ing that the density matrix be of rank 2, give rise to very 
peculiar dynamics of the quantum correlations in the time- 
evolved states p^g under the corresponding ^-type channels. 
Namely, defining the threshold time yt^ = - 1 In c, we find an- 
alytically that every measure Q considered in this paper takes 
the form 



Qipf^'^ - ^ ^^^'^ 






(3) 



This entails that all studied measures of discord rcmsiin frozen 
to their initial value until, at ^ = t^, they suddenly start de- 
caying exponentially, as depicted in Fig. [T] The functional 
dependence fa for each measure Q is reported in Table |l] 

The presented results, although tied to the particular 
choices of channels and initial states, apply to such a variety 
of measures of discord that the level of mere coincidence can 
be safely considered overcome. We now seek for a physical 
explanation of our findings. Most of our employed measures, 
in the chosen settings, can be interpreted as minimal distances 



from the set of CQ states. It turns out that, for nondissipa- 
tive evolutions with initial condition as in Eq. ([2]), all our dis- 
tance functions are optimized by the same time-dependent CQ 
state (explicitly given in |[24l[26l[66l). This certainly explains 
why all such measures have a sudden change at the same time 
t = t^, but does not explain why all of them are in fact con- 
stant for earlier times. For instance, without the purity ad- 
justment, the geometric discord Dg 144 1 decreases rather than 
being frozen when t < t* i66ll68l. The adjusted Dq instead 
freezes like the other measures. As remarked earlier, the pu- 
rity adjustment is necessary to correct for some flaws of Dg 
||46l|5P|. This suggests that the occurrence of freezing is, at 
a somewhat empirical level, a stronger feature which necessi- 
tates to be present in bona fide quantifiers of discord. 

A deeper insight into this phenomenon can then be achieved 
by looking to the geometric representation of BD states, sim- 
ilarly to what was done in f30] and in [68] for the quantum 
discord D and for the (non-adjusted) geometric discord Dq, 
respectively. In the space (ci,C2,C3), one can draw the sur- 
faces along which each given measure Q is constant, and then 
superimpose the trajectories corresponding, say, to the evolu- 



Measure Q 



fa(s) 



Quantum discord V 

One-way quantum deficit A^ 

Relative entropy of discord Dr 

Adjusted geometric discord Dg 

Trace-distance discord Di 

Negativity of quantumness Q^ 

Bures-distance discord Db 

Fidelity-based measure Dp 

Local quantum uncertainty "U 



2^2/(1 + s^) 



(1 + V2)[i-(vr 
1- vi-52 



'■ + iy] 



TABLE I: Summary of the expressions for the measures Q of quan- 
tum correlations in the BD states with initial condition c^^^(O) evolv- 
ing under local nondissipative channels, as evaluated in Eq. ([3]). 



tion of c^^\t) under local ^-channels, in order to visualize the 
freezing. Let us focus, without loss of generality, on phase 
flip channels (k = 3). For ease of graphical display, instead 
of drawing three-dimensional surfaces, we find it more infor- 
mative to slice the tetrahedron of BD states at, say, a constant 
value of C3, and draw the contours of constant Q in the re- 
sulting projected plane spanned by (ci,C2). This analysis is 
reported in Fig. [2] for several types of discord measures. It is 
evident that the various quantifiers exhibit their quite diff'erent 
nature even within the restricted set of BD states. In partic- 
ular, around the classical states at the centre (ci, C2) = (0, 0), 
contours of constant quantum correlations are ellipses for en- 
tropic (D, A~^ , Dr) and fidelity -based (DB,Df,^) measures, 
while they are circles for the Hilbert- Schmidt based Dq, and 
squares for the trace-distance based Di. Similarly, closer to 
the peripheral boundaries of the physically allowed region, 
the contour lines have diff'erent topologies across the various 
measures. Some measures like Di are constant for wider re- 
gions than the others, for instance. Here we are interested in 
the overlap between the contour lines of all (3's. We find that 
the only possible straight lines, in the BD state space, which 
keep all the considered measures simultaneously constant, are 
those along the diagonals of the projected rectangles in the 
space (ci, C2), stopping sufficiently far away from the central 
core (see Fig.|2]) |62|. These segments describe only and pre- 
cisely the evolutions of BD states with initial conditions c^^\0) 
[Eq. ([2])] for local phase flip channels — and similarly for the 
other ^-channels by selecting the corresponding planes (c/, Cj) 
and conditions ^^^(0) — up to the time t*. This yields a general 
geometric interpretation to the phenomenon of frozen discord 
universally observed in all valid measures at once. 

We wish to remark that, in the presence of suitable non- 
Markovian channels which can be described by a master equa- 
tion with a memory kernel (as in the case of pure dephasing 
or decoherence under classical random external fields) 125] - 
[28l [661 . all the phenomena previously observed for the en- 
tropic discord D extend automatically to all the measures con- 
sidered in this paper. Indeed, in the models of i25H28ll66l , the 
evolution of BD states can be formally written as in Eq. ([T]), 
but with 2yt replaced by a more general A(0, which could 
be decreasing over some time intervals. This can give rise 
to a dynamics with multiple intervals of constant discord 
(2611281 [66 1, or discord frozen forever |25| depending on the 
initial conditions. By our analysis, we conclude that those fea- 
tures, which might be observable e.g. in the dynamics of im- 
purity atoms in Bose-Einstein condensates 1251 169|, are uni- 
versal too and detectable by any suitable discord measure Q. 

In this paper we established the general status of an intrigu- 
ing aspect of quantum correlations other than entanglement: 
the freezing in certain dynamical conditions. Originally re- 
vealed for entropic quantifiers 1241 . we showed that this fea- 
ture is common to all bona fide measures of discord, and we 
provided a geometric interpretation thereof. The occurrence 
of freezing in specific nondissipative evolutions can be pro- 
posed as a 'sanity check' to validate novel discord quantifiers. 
From an operative perspective, we showed that quantumness 
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FIG. 2: (Color online). Contours of constant Q in the plane (ci, C2) 
for physical BD states with C3 = ^, where Q represents (a) T> = 
A^ = Dr, (b) Dg, (c) Di = Qn, (d) Db. The diagonals (traced from 
the vertices to the centre) represent the evolution of BD states with 
initial condition c^^^(O) [Eq.|2J under local phase flip channels. For 
t < f*" (continuous traits) all the measures Q are constant; they then 
decay to zero for later times (dashed traits). The color legend for the 
(3's in the contour plots is: ^^^H 1. 



of correlations, in all its manifold manifestations, can be sus- 
tained with no loss in suitable noisy settings. Our predictions 
are amenable to experimental verification with current tech- 
nology, e.g. using photons |70 1 or nuclear magnetic resonance 
techniques |71, 72 1. We expect our study to stimulate novel 
endeavors in the comprehension and exploitation of genuinely 
quantum eff'ects in open systems. 
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SUPPLEMENTAL MATERIAL 

In this supplemental material we calculate the fidelity- 
based measure of discord for arbitrary Bell diagonal states, 
and we prove that, in the state space of Bell diagonal states, 
the only straight-line dynamical trajectories, under which all 
studied measures of discord freeze, are those corresponding to 
the initial conditions specified in the main text. 

Fidelity-based measure of discord 

A two-qubit state p can be written in the Bloch representa- 
tion: 

-.3 3 3 

(A.l) 
where {crf'^} denote the Pauli matrices. With such a state p 
we associate the triple {x,y, T} Ii44il . 

Theorem 1. For any Bell diagonal (BD) state p, 
F(p,crQ) > F(p,cr), where cr is the state with triple 
{x, y, T}, and ctq the state with triple {0, 0, T}. 

Proof. For any state cr, there exists another state with asso- 
ciated triple {-X, -y, T} which we denote cr_ . It can be verified 
that for BD p, the matrices A/pa ^|p and ^fpcr- ^ have the 
same characteristic polynomial, and so have the same eigen- 
values. From the expression for fidelity 



For any state in this form, a second state ;^o = p'WiM^i\^ ® 
p'^ + (1 - p')Mfik^i\^ ®p2 can be derived using the identity 



F(p, cr) = tr 



y^o-^ 



(A.2) 



It can be seen that the equality of the eigenvalues is sufficient 
to prove the identity 



F(p,cr) = F(p,cr_) 
Fidelity obeys the concavity property 



(A.3) 



F(p,yUO-i -h (1 -yU)0-2) > yUF(p, O-i) -h (1 -yU)F(p, 0-2),yU G [0, 1] 

(A.4) 



Setting \i - \,cr\ = o-, 0-2 = o--, using Eq. (A.3) and noting 



that ^(cr -h (T_) = (To, we recover F(p,o-o) > F(p,o-) and so 
prove the theorem. n 

Equipped with this, it is now possible to prove the main 
result. 



Theorem 2. For any BD state p, the classical-quantum 
(CQ) state X which maximizes F(p,x) is also BD 

Proof. Any CQ two-qubit state will be of the form x = 
/^l^iX^il^ ^pf + (1 - /?)l^2X^2l^ ^P^ where {J^Ai)^, 1^2)^} is 
an orthonormal basis for qubit A. Such a CQ state will have 
associated triple {(2p - l)e, s+, e^_} |44|, where 



et = <<Aik#i>, (A.5) 

s±j = tr [(ppi ± (1 - p)p2)cri] . (A.6) 



"■4 

1^2)'^ = 1^2)^ 

pf = i(i^ + K-a-/^)Tf) 



(A.7a) 

(A.7b) 

(A.7c) 

(A.7d) 

(A.7e) 



where ti, T2 are the traceless parts of pi,p2. This state is man- 
ifestly CQ, and it can be easily verified that it will have as- 
sociated triple {0, 0, T}. From Theorem 1. we can see that it 
suffices to restrict ourselves to CQ states of this form. 

Temporarily, we relax the restriction that ^is of unit length, 
and consider the set of states where \\e\\ < 1. This allows 
us to repeat the previous trick, this time between x with e = 
(^1, ^2, ^3), si = (si, S2, 0) din&x' with e' - (^1, ^2, -^3), s'_ - 
s"- . As before, from the comparison of characteristic polyno- 
mials, F{p,x) - ^(P^X')- A similar result holds for ^1 = 
and ^2 = and, by switching the vectors we consider, for any 
et = 0. From this we can see that if any si = then for maxi- 
mum fidelity, ei = 0, and vice versa. We also note that this set 
of states is convex, and so due to the concavity of fidelity, any 
local maximum will be a global maximum. 

As an ansatz, we now consider the states where 
ej = eSij, Sj = sSij where / sets the non-zero vector 
element. From the previous result we can see that maximiza- 
tion only needs to be performed over ei and si as fidelity 
can only decrease under any variation in any single other 
element. Furthermore, et and st appear only as a product 
efSi in the density matrix, never on their own. This means 
that maximizing over both is equivalent to setting et = 1 and 
maximizing only over si, thus allowing us to reimpose the 
restriction that \\e\\ = 1 and returning to the CQ states. It is 
now only necessary to maximize over a single parameter, 
and any local maximum in this parameter will be the global 
maximum for CQ states. The remaining states are the BD 
CQ states, and so finding the maximum among these states 
proves the theorem. 

Maximizing over the single remaining parameter we obtain 
the result for BD states with eigenvalues a,fi, 7, S: 



where 



Df(p) = min{quq2,q3} (A.8) 

^1 = 1 - V2 ( V^ -h ^/^) (A.9a) 

q2 = l- V2 ( V^ -h ^|^) (A.9b) 

q3 = l- V2 ( V^ -h Vr^) (A.9c) 



This is identical to the expression for local quantum uncer- 
tainty [,18,1 of BD states, giving the identity 

Df(p) = 'U(p) (A. 10) 

forp an arbitrary BD two-qubit state. n 

Universal freezing trajectories 

Theorem 3. Universal freezing for decay in the k-type flip 
channel is observed in BD states only for initial conditions of 
the form 



Jk) 



cTXto) = ±e 



-2yto 



, cy(to) 



+c['\to)e-^^', with \c['\to)\ ^ c 
(A. 11) 



Proof We consider a BD state in the formp^B = ^C^ab + 
X^zzi ^/<^ ^ erf). For universal freezing, all bona fide mea- 
sures must have a constant value up until the time t* . In par- 
ticular, any two of them must satisfy this requirement. In the 
following, we will pick the negativity of quantumness and the 
adjusted geometric discord and find that the only straight-line 
freezing trajectory common to both is the one reported in the 
claim. Recall the expression for trace-distance discord (Di) 
alias negativity of quantumness (Qn), Di = Qn = ^2, where 
{gt} are simply {|c/|} ordered such that ^1 > ^2 > ^3. From this 
it can immediately be seen all initial BD states p satisfying 



^2(^0) 



\cf\to)\ ^ c 



(A.12) 



will exhibit freezing for these measures. 

The adjusted geometric discord Dq for the same states has 
the following expression for time t < t^: 



Dg = 2{c^ + ^l)l{\ + ^i + c^ + e^) 



(A.13) 



Since ^/(O = Q(to)e~^^^ for / = 1,3, time evolution will 
follow a straight line £ in the f 1 - ^3 plane. Such a line can be 
defined parametrically as: 



^^(s) = ms -\- a 



(A. 14a) 
(A. 14b) 



In order for Dq to be constant over this time interval, £ must 
follow a contour line in this plane. Since contour lines are 
perpendicular to the gradient, it must therefore be: 



VDg.V^ = 



(A. 15) 



Calculating this and then substituting in the values for f 1 and 
^3 from Eq. ( |A.14| ) gives: 



a^s -\- am (s^-l]-\-s(c^- m^\ = 



(A. 16) 



In order for ^ to be a straight line, m and a must be independent 
of ^. It can readily be verified by solving for a that the only 
solution for which this is the case is a = 0, m = ±c. 



Points on these lines correspond exactly to Eq. ( |A.11| ), 
proving the theorem. n 



